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Abstract 

A  three-dimensional  mathematical  model  for  PEM  fuel  cells  operated  on  reformate  is  developed  based  on  our  previous  established  fuel  cell 
model  [Int.  J.  Transport  Phenomena  3  (2001)  177],  by  incorporating  the  adsorption  and  oxidation  kinetics  of  CO  on  platinum  surface  proposed 
by  Springer  et  al.  [Proceedings  of  the  Electrochemical  Society,  Montreal,  Canada,  1997;  J.  Electrochem.  Soc.  148  (2001)  All].  This  model 
is  capable  of  studying  the  effect  of  CO  poisoning  as  well  as  the  hydrogen  dilution  effect  by  inert  gases.  The  adsorption  and  oxidation  kinetics 
of  CO  on  a  platinum  surface  are  incorporated  in  the  source  terms  of  the  species  equations;  thus,  the  basic  form  of  the  mathematical  equations 
are  the  same  as  those  used  for  PEM  fuel  cells  operated  on  pure  hydrogen.  With  this  model,  we  can  obtain  detailed  information  on  the  CO 
poisoning  and  variation  of  CO  and  hydrogen  concentrations  inside  the  anode.  The  results  from  this  3D  model  reveal  many  new  phenomena 
that  cannot  be  obtained  from  previous  ID  or  2D  models.  Results  of  the  effects  of  various  operating  and  design  parameters,  such  as  anode 
flow  rate,  gas  diffuser  porosity,  gas  diffuser  thickness,  and  the  width  of  the  collector  plate  shoulder,  are  also  presented.  The  modeling  results 
demonstrate  the  value  of  this  model  as  a  design  and  optimization  tool  for  the  anode  of  PEM  fuel  cells  operating  on  reformate. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  hydrogen  PEM  fuel  cells,  when  pure  hydrogen  is  used 
as  the  fuel,  usually  very  low  polarization  occurs  at  the  anode 
due  to  the  fast  hydrogen  oxidation.  Thus  the  overpotential 
at  the  anode  is  often  neglected  in  fuel  cell  models  [4-6]. 
If  reformate  is  used  as  the  fuel  for  PEM  fuel  cell,  the  di¬ 
lution  effect  of  the  inert  gases  at  the  anode  can  no  longer 
be  neglected.  More  importantly,  traces  of  carbon  monoxide 
(CO)  are  inevitably  present  in  the  reformed  gases  and  it  can 
cause  substantial  poisoning  of  the  catalyst  and  lead  to  dras¬ 
tic  degradation  of  fuel  cell  performance  [7].  CO  affects  the 
reactivity  of  the  catalyst  mainly  by  preventing  H2  adsorption 
as  the  CO  molecules  adsorb  strongly  on  the  active  sites  of 
the  catalyst  surface  [8].  This  site  exclusion  behavior  of  CO, 
combined  with  the  dilution  effect  due  to  decreased  hydrogen 
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partial  pressure  in  the  anode  channel,  can  cause  significant 
cell  voltage  losses.  These  losses,  when  added  to  the  losses  at 
the  cathode,  can  reduce  the  fuel  cell  efficiency  to  an  unus¬ 
able  level.  Ways  to  reduce  the  CO  poisoning  effects  include 
internal  cleanup  methods,  such  as  oxygen  bleeding  and  other 
oxidation  means  [9],  using  catalyst  material  resistant  to  the 
poisoning  [10],  and  adopting  operation  techniques  that  can 
improve  CO  tolerance  [11].  Other  research  includes  studies 
on  the  mechanism  of  CO  poisoning  [12-16]  and  the  CO  tol¬ 
erance  [17,18].  The  performance  data  from  PEMFC  operated 
on  reformed  gas  can  be  found  in  Refs.  [20,21].  The  kinetic 
model  presented  by  Springer  et  al.  [2,3]  succeeded  in  explain¬ 
ing  many  features  of  the  CO  adsorption  and  reaction  during 
the  PEMFC  operation.  To  account  for  the  changing  of  CO 
inside  the  fuel  cell  chamber,  Zhang  et  al.  [21]  presented  a  CO 
inventory  model,  which  used  the  basic  mass  balance  between 
the  inlet  and  outlet  to  calculate  the  consumption  of  CO  in¬ 
side  the  fuel  cell.  Though,  this  model  calculated  the  CO  mole 
fractions  at  the  anode  feed  and  the  anode  outlet,  it  neglected 
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Nomenclature 

(a/QCf)a  reference  exchange  current  density  times  area 
at  anode  (A/m3) 

(a/r0ef)c  reference  exchange  current  density  times  area 
at  cathode  (A/m3) 

Ag  effective  activate  surface  area  per  unit  volume 
of  the  catalyst  layer  (m2/m3) 
bfc  back-to-forward  CO  adsorption  ratio  (atm) 

Z?fh  back-to-forward  H2  adsorption  ratio  (atm) 
c  mole  concentration  of  the  gas  mixture 

(mol/m3) 

cH+  proton  mole  concentration  (mol/m3) 

cp  specific  heat  at  constant  pressure  (J/kg  K) 

Dk  binary  diffusion  coefficient  of  the  Cth  compo¬ 
nent  (m2/s) 

E  fuel  cell  voltage  (V) 

Eq  open  circuit  potential  (V) 

F  Faraday  constant,  96  487  (C/mol) 

i  current  density  (A/m2) 

ico  local  kinetic  current  density  from  CO  (A/m2) 

/h  local  kinetic  current  density  from  hydrogen 
(A/m2) 

j  transfer  current  density  (A/m3) 

k  thermal  conductivity  (W/m  K) 

kQC  pre-exponential  of  CO  Tafel  electro-oxidation 

rate  (A/m2) 

keh  pre-exponential  of  //a d  Tafel  electro-oxidation 

rate  (A/m2) 

kfc  electrode  forward  H2  adsorption  rate  times  2F 
(A/cm2  atm) 

kfh  electrode  forward  CO  adsorption  rate  times  2F 
(A/m2  atm) 

kh  hydraulic  permeability  of  the  porous  layer  (m2 ) 

ky  electrokinetic  permeability  (m2) 

lct  catalyst  layer  thickness  (m) 

P  pressure  (Pa) 

Q  heat  generation  rate  of  the  energy  equation 

r(2)  correction  factor  in  the  momentum  equations 
for  the  porous  media 

R  universal  gas  constant  (8.314  J/mol  K) 

Sk  source  term  of  the  species  transport  equation 

T  temperature  (K) 

V  velocity  vector  (m/s) 

X  mole  fraction 

Greek  letters 

0  electrode  overpotential  (V) 

s  porosity 

fj  viscosity  (kg/m  s) 

cp  membrane  phase  potential  (V) 

p  density  (kg/m3) 

a  ionic  conductivity  (£2_1  m-1) 


$co  fractional  CO  coverage  of  Pt  surface 
$h  fractional  coverage  of  Pt  surface,  hydrogen 

Subscripts 
a  anode 

c  cathode,  CO 

ct  catalyst  layer 

eff  effective 

h  hydrogen 

i  inlet 

k  k’th  component 

m  membrane 

p  porous  media 

Superscripts 
a  anode,  anodic 

c  cathode,  cathodal 

ref  reference 


the  variation  of  CO  concentration  along  the  anode  flow  direc¬ 
tion  as  well  as  in  across-the-channel  direction.  In  addition, 
the  changing  of  hydrogen  concentration  across  the  fuel  cell 
and  along  the  flow  direction  was  not  taken  into  account. 

In  this  paper,  the  kinetic  equations  for  the  CO  adsorp¬ 
tion  and  oxidation  behavior  provided  by  Springer  et  al.  [2] 
is  incorporated  into  our  previously  established  3D  model  to 
study  the  CO  poisoning  effect  in  a  fuel  cell.  In  addition,  the 
dilution  effect  by  the  inert  gases  is  also  included.  With  this 
model,  we  can  obtain  a  detailed  variation  of  CO  and  hydrogen 
concentrations  inside  the  anode.  Unlike  the  previous  models 
[2,3,21],  the  variation  of  the  chemical  reaction  inside  the  fuel 
cell  is  fully  accounted.  The  modeling  results  also  compared 
very  well  with  the  experimental  data  by  Divisek  et  al.  [19]. 


2.  Mathematical  model 

Fig.  1  is  the  schematic  of  the  model  geometry  of  a  single 
fuel  cell.  The  details  of  the  general  3D  model  can  be  found 
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Fig.  1.  The  schematic  of  the  model  geometry  of  a  single  fuel  cell. 
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in  Zhou  and  Liu  [1].  A  summary  of  the  model  equations  is 
provided  below: 

Continuity  equation: 

V(pV)  =  0  (1) 

Momentum  equation: 


The  heat  source  Q  is  given  by, 

f  ljn 

J  '^a  anode  catalyst  layer 


Q={ 


•2 


<Xi 


membrane 
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^  ^C°ci  1  cathode  catalyst  layer 
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pV(VY)  =  -  VP  +  < 
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r®/iV2Y 


r(2VV2V 


£/Z 

kfV 


S/JL 
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+  £Cjj+ 


kp 


gas  channel 
gas  diffuser 

catalyst  layer  and  membrane 


(10) 


Species  equation: 


The  membrane  phase  potential  satisfies: 


pV(VXfc)  =  spDt  nV2Xk 


+ 


0  gas  channel  and  gas  diffuser 
spSk  catalyst  layer 


Energy  equation: 
pcpV(V7yO  =  ^effV27^ 


Jc 


V(crV?7)  = 


cathode  catalyst  layer 

membrane 

anode  catalyst  layer 


Then  the  cell  voltage  is  determined  by: 


E  =  Eo\r]c\  -  1 77a I  - 


m 


CTmdj  J  ct.aCTctdy 


(ID 


+ 


0  gas  channel  and  gas  diffuser 
Q  catalyst  layer  and  membrane 


The  correction  parameters  A2\  Deff  are  defined  as  [23]: 


r(2)  =  2.25 


(5) 


Dk  gas  channel 

Deff  -  ^  j  r 

DkS  porous  media 


(6) 


^eff  — 


"gas 


gas  channel 


~2kP  +  g/(2fep+A:gas)+( l  — g)/3A:p  porous  media 


(7) 


At  the  cathode,  the  mass  generation  source  term  S&,  for  oxy¬ 
gen,  water,  and  protons  are  jJ(4Ec ),  —jc/(2Fc),  and  jJ(Fc ), 
respectively.  At  the  anode,  the  source  terms  for  hydrogen 
molecules  and  protons  are  — jJ(2Fc ),  and  jJ(Fc ),  respec¬ 
tively.  The  volumetric  current  generation  rates,  ja  and  jc,  are 
calculated  from  Butler- Volmer  expression: 


(12) 


where  Eq  is  the  open  circuit  potential.  In  this  paper,  the  cor¬ 
relation  obtained  by  Gurau  et  al.  [22]  from  the  experimental 
data  of  Parthasarathy  et  al.  [23]  is  used  for  the  specific  oper¬ 
ation  conditions  listed  in  Table  1 : 


E0  =  0.00257  +  0.2329  (13) 

Note  that  generally  Eq  depends  on  pressure  as  well  as  temper¬ 
ature,  but  in  this  study,  all  the  cases  have  the  same  operating 
pressure.  In  Eq.  (12),  the  values  of  the  integrals  depends  on  x 
and  z  since  current  density  is  a  function  of  x  and  z,  in  addition 
to  y.  In  determining  the  cell  voltage,  the  average  value  of  each 
integral  is  used,  thus  inexplicitly  taking  into  account  of  the 
current  re-distributions  in  GDL  on  the  x-z  plane. 

In  the  presence  of  CO  at  the  anode,  part  of  the  catalyst 
surface  is  covered  by  CO  that  causes  site  exclusion  and  elec¬ 
tron  capture  through  dipole  interactions.  The  value  of  (tfJoef)a 
used  in  Eq.  (8)  is  then  reduced.  Spring  et  al.  [2,3]  suggested 
that,  in  the  presence  of  CO,  the  interfacial  kinetics  of  PEMFC 
using  reformate  could  be  represented  by  four  processes: 


ja=K)e,)a(— y/2 

V  Ct2ref  / 


X 


exp 


<E_ 

RT 


?7a 


exp 


<E_ 

RT 


?7a 


Jc  = 


\  cC>2ref  / 


X 


exp 


<1 

RT 


0c 


exp 


<E_ 

RT 


0c 


(8) 


(9) 


(1)  CO  adsorption: 

^fc 

CO  +  Pt  ^  (Pt  -  CO) 

bfckfc 

(2)  Dissociative  chemisorption  of  hydrogen: 

^fh 

H2  +  2Pt  <=►  2(Pt  —  H) 

^fh^fh 

(3)  Electrochemical  oxidation  of  hydrogen: 
(Pt  -  H)  H+  +  e~  +  Pt 


(14) 


(15) 


(16) 
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Table  1 

Base-case  electrode  parameters  and  propertiesa 

Gas  channel  length,  8.0  x  10-2  m 
Gas  channel  height,  8.0  x  10 * *  4  m 
Gas  diffuser  thickness,  2.0  x  10  4  m 
Catalyst  layer  thickness,  0.3  x  10  4 *  m 
Gas  channel  width,  10  x  10-4  m 
Collector  plate  shoulder  width,  10  x  10  4  m 
Membrane  thickness,  2.3  x  10-4  m 
Gas  diffuser  porosity,  s  =  0.4 
Membrane  porosity,  sm  =  0.28 

Volume  fraction  of  membrane  in  catalyst  layer,  £mc  =0.5 
Permeability  to  air  of  the  gas  diffuser,  kp  =  1.76  x  10  11  m2 * * * * * 
Hydraulic  permeability  of  membrane,  km  =  1.58  x  10  18  m2 
Reference  exchange  current  density  times  area  at  the  cathode,  (ailf  )c 
=  5  x  102  A/m3 

Oxygen  reference  concentration,  3.39  mol/m3 

Hydrogen  reference  concentration,  56.4  mol/m3 

Air  inlet  pressure,  3  atm 

Air  inlet  velocity,  3  m/s 

Air  inlet  temperature,  82  °C 

Oxygen/nitrogen  ratio  in  dry  air  at  inlet,  0.21/0.79 

Cathode  Reynolds  number  =  347.68 

Anode  inlet  pressure,  3  atm 

Anode  inlet  velocity,  3  m/s 

Anode  inlet  temperature,  82  °C 

Anode  Reynolds  number  =  175.78 

Anode  inlet  composition  on  dry  base:  60%  H2,  40%  N2/CO2,  50  ppm  CO 

Ionic  conductivity  of  membrane,  crm  =  10  £2~[  m-1 

Reference  exchange  current  density  times  area  at  the  cathode,  (a?Qf)a= 

5  x  108 * *  A/m3 

Effective  area  per  unit  volume  of  the  catalyst  layer  Ae  =  50  000  m2/m3 
Back-to-forward  CO  adsorption  ratio,  bfc  =  2  x  10  5  atmb 
Back-to-forward  H2  adsorption  ratio,  b^  =  0.5  atmb 
Electrode  forward  H2  adsorption  rate  times  2F,  k\c  =  2.4  x  105  A/(m2  atm)b 
Electrode  forward  CO  adsorption  rate  times  2F,  kfh  =  40  000  A/(m2  atm)b 
Pre-exponential  CO  Tafel  electro-oxidation  rate,  kec  =  10“2  A/m2b 
Pre-exponential  Ha(j  Tafel  electro-oxidation  rate,  ke  1,  =  5000  A/m2b 

a  Most  data  from  Gurau  et  al.  [22]  and  Zhou  and  Liu  [1]. 

b  From  or  adapted  from  Springer  et  al.  [3]. 

(4)  Electrochemical  oxidation  of  CO  to  CO2: 

H20  +  (Pt  -  CO)  A  Pt  +  C02  +  2H+  +  2e~  (17) 

According  to  Springer  et  al.  [2,3],  the  first  two  processes 

occur  simultaneously  to  compete  for  the  vacant  catalyst  sur¬ 

face.  The  third  equation  represents  the  current  generation  by 

the  anodic  oxidation  of  hydrogen  and  the  fourth  corresponds 

to  the  electrochemical  oxidation  of  CO  to  CO2,  which  reaches 

significant  rate  only  at  high  anode  overpotential.  The  poi¬ 

soning  effect  of  CO  is  mainly  due  to  the  chemical  reaction 
represented  by  Eq.  (14).  Due  to  this  process,  the  part  of  the 

Pt  surface  that  adsorbs  CO  is  unavailable  for  the  chemical 
reaction  represented  by  Eq.  (15)  and  subsequent  process  by 
Eq.  (16).  At  large  anode  overpotentials,  the  chemical  reaction 
represented  by  Eq.  (17)  will  reach  a  significant  level.  Thus 

CO  will  be  consumed  and  part  of  the  CO-adsorbed  Pt  surface 

will  then  be  freed. 

Under  steady-state  assumptions,  Springer  et  al.  [2]  solved 

the  CO  and  hydrogen  adsorption/desorption  equilibrium 

equations  and  obtained  the  local  generated  ionic  current 


densities  as: 

/h  =  2fceiAsinh(Aj  (18) 

ico  =  2kec6Co  e‘h/hc  (19) 

where  b h  =  2RT/F  and  bc  =  2RT/F. 

Then  the  transfer  current ja  within  the  anode  catalyst  layer 
becomes: 

Ja  =  tfaO'h  +  z'co)  (20) 

At  the  anode,  the  consumptions  of  hydrogen  and  CO  are: 

-Ae6i  -Ae/co  nu 

2  Fc  ’  2  Fc  (  j 

where  Ae  is  the  effective  active  catalyst  surface  area  per  unit 
volume  of  the  catalyst  layer. 

As  the  current  is  mainly  from  the  anode  to  the  cathode, 
the  components  of  the  current  along  the  x  and  z  axes  (Fig.  1) 


Fig.  2.  Modeling  results  with  varied  CO  content,  (a)  current  density  vs. 
overpotential  (b)  polarization  curve.  Case  conditions:  anode  flow  is  fully 
saturated  reformed  gases,  with  60%  hydrogen,  40%  N2/CO2  on  dry  base, 
flow  velocity  =  1  m/s,  temperature  =  82  °C,  pressure  =  3  atm. 
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are  neglected,  and  we  have, 
di 

—  =  0aOh  +  *co)  at  anode  (22) 

dy 

Because  i  =  cr(drj/dy)  and  no  protons  pass  through  the  inter¬ 
face  between  the  GDL  and  the  catalyst  layer,  thus: 

dri 

=  0  (23) 

dy 

Assuming  constant  overpotential  along  the  interface  between 
the  catalyst  layer  and  the  membrane,  we  have: 

0  =  da  (24) 

For  the  above  mathematical  model,  inlet  flow  conditions  are 
prescribed  at  the  gas  channel  inlets.  The  boundary  conditions 
for  the  species  transport  equations  include  the  Neumann  con¬ 
ditions  prescribed  along  the  gas  channel  walls  as  well  as  along 
the  interface  between  the  catalyst  layer  and  the  membrane. 
For  the  energy  equation,  constant  temperatures  along  the  gas 


channel  outer  walls  are  used  to  simulate  bath  cooling.  Finally, 
between  the  different  layers  within  the  fuel  cell,  continuity  is 
assumed  at  the  interfaces.  With  a  computer  program  devel¬ 
oped  in-house  [1],  the  governing  equations  in  the  different 
regions  of  the  fuel  cell  are  solved  simultaneously  to  ensure 
the  coupling  of  the  flow,  species,  overpotential  and  current 
density. 

3.  Computational  results  and  discussions 

A  series  of  computation  simulations  were  performed  using 
the  3D  fuel  cell  model.  The  base  case  data  and  parameter 
values  are  listed  in  Table  1 . 

Fig.  2(a)  shows  the  modeling  results  of  anode  overpoten¬ 
tial  changes  with  different  current  densities  at  different  inlet 
CO  concentrations,  and  Fig.  2(b)  shows  the  corresponding 
polarization  curves.  From  Fig.  2(a),  it  is  seen  that  there 
is  substantial  overpotential  at  the  anode  when  even  a  very 


0.472 

0.446 

0.419 

0.393 

0.367 

0.340 

0.314 

0.288 


(b) 


(b) 


flow  direction 


Fig.  3.  Hydrogen  mole  fractions  along  (a)  flow  direction  (b)  across  flow 
direction.  Case  conditions:  anode  flow  is  fully  saturated  reformed  gases,  with 
60%  hydrogen,  40%  N2/CO2  on  dry  base,  flow  velocity  =  1  m/s,  temperature 
=  82  °C,  pressure  =  3  atm,  average  current  density  =  0.64  A/cm2,  anode 
stoichiometric  ratio  =  3.39. 


Fig.  4.  CO  mole  fractions  along  (a)  flow  direction  (b)  across  flow  direction. 
Case  conditions:  anode  flow  is  fully  saturated  reformed  gases,  with  60%  H2, 
40%  N2/CO2,  50 ppm  CO  on  dry  base,  flow  velocity  =  1  m/s,  temperature 
=  82  °C,  pressure  =  3  atm,  average  current  density  =  0.64  A/cm2,  anode 
stoichiometric  ratio  =  3.39. 
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small  amount  of  CO  is  present.  The  larger  the  CO  content 
is,  the  greater  is  the  anode  overpotential.  As  a  result,  the 
fuel  cell  performance  degrades  drastically  as  shown  by  the 
polarization  curves  in  Fig.  2(b).  From  Fig.  2(a),  it  is  also 
seen  that  at  low  current  densities  the  CO  poison  effect  is 
not  as  significant.  This  CO-tolerance  exists  because  only  a 
small  portion  of  the  catalyst  surface  is  needed  to  generate 
the  currents.  At  such  low  current  densities,  the  reduction  of 
the  active  area  of  the  catalyst  layer  does  not  have  significant 
effect  on  the  rate  of  hydrogen  oxidation  yet,  due  to  the 
very  fast  hydrogen  oxidation  kinetics.  As  current  density 
increases,  more  active  area  of  the  catalyst  surface  is  needed. 
Then  the  reduction  of  the  active  area  due  to  CO-adsorption 
by  the  Pt  surface  has  significant  effect  on  the  rate  of  hydro¬ 
gen  oxidation.  To  compensate  for  the  CO-covered  catalyst 
surface,  larger  activation  loss  at  the  anode  incurs  as  shown 
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by  Fig.  2(a).  This  rapid  increase  of  overpotential  continues 
until  a  certain  level  when  the  rate  of  this  increase  is  reduced. 
At  this  point,  the  CO  oxidization  reaches  a  significant  level 
and  the  consumption  of  CO  releases  some  CO-covered 
Pt.  Eventually  the  curves  in  Fig.  2(b)  become  parallel, 
indicating  that  the  absorption  and  oxidation  of  CO  reach 
equilibrium. 

Fig.  3(a)  and  (b)  are  the  typical  concentration  contours  of 
hydrogen  mole  fraction  along  the  flow  direction  and  across 
the  fuel  cell.  Fig.  4(a)  and  (b)  are  the  corresponding  CO 
concentration  contours.  Along  the  flow  direction,  as  CO  and 
hydrogen  is  oxidized,  CO  and  hydrogen  concentrations  are 
reduced.  A  comparison  of  the  local  current  densities  of  % 
and  ico  shows  that  4  is  far  larger  than  /co-  This  indicates 
that  CO  oxidation  affects  the  PEMFC  performance  through 
the  release  of  Pt  surface  that  is  previously  covered  by  CO, 
while  the  contribution  of  the  current  density  generated  by  the 
CO  oxidation  is  negligible. 


Fig.  5.  Local  current  density  along  (a)  flow  direction  (b)  across  flow  direc¬ 
tion.  Case  conditions:  anode  flow  is  fully  saturated  reformed  gases,  with 
60%  H2,  40%  N2/CO2,  50  ppm  CO  on  dry  base,  flow  velocity  =  1  m/s,  tem¬ 
perature  =  82  °C,  pressure  =  3  atm,  average  current  density  =  0.64  A/cm2, 
anode  stoichiometric  ratio  =  3.39. 


Fig.  6.  Current  density  vs.  overpotential  and  polarization  curves  when  the 
current  densities  are  averaged  at  the  fuel  cell  inlet,  outlet  and  whole  fuel 
cell  chamber.  Case  conditions:  anode  flow  is  fully  saturated  reformed  gases, 
with  60%  H2,  40%  N2/CO2,  50  ppm  CO  on  dry  base,  flow  velocity  =  1  m/s, 
temperature  =  82  °C,  pressure  =  3  atm. 
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A  new  phenomenon  from  the  3D  modeling  effort  is  that, 
due  to  the  reduction  of  CO,  near  the  exit  of  the  fuel  cell, 
the  chemical  reaction  rate  is  actually  greater  than  that  near 
the  inlet.  This  is  obviously  not  the  case  for  a  fuel  cell  oper¬ 
ated  on  pure  hydrogen.  This  phenomenon  is  shown  by  the 
local  current  density  distribution  in  Fig.  5(a)  and  (b).  Near 
the  exit,  though  hydrogen  concentration  is  reduced  due  to 
reaction  consumption,  the  fast  oxidation  of  hydrogen  makes 
the  influence  of  reduced  hydrogen  concentration  negligible. 
On  the  other  hand,  the  lower  CO  concentration  down  the  flow 
channel  due  to  CO  oxidation  reduces  the  poisoning  effect  and 
leads  to  higher  current  density.  To  see  this  reduction  effect 
of  CO  poisoning  more  clearly,  the  relationships  between  cur¬ 
rent  density  and  anode  overpotentials  at  the  inlet  and  outlet 
are  plotted  in  Fig.  6(a)  and  compared  with  the  average  value 
for  the  cell.  The  corresponding  polarization  curves  averaged 
at  the  inlet  and  the  outlet  as  well  as  for  the  whole  cell  are 
plotted  in  Fig.  6(b). 


Fig.  7.  Current  density  vs.  overpotential  and  power  output  vs.  current  density 
when  the  current  densities  are  averaged  at  regions  over  the  shoulder,  over 
the  channel  and  the  whole  fuel  cell  chamber.  Case  conditions:  anode  flow  is 
fully  saturated  reformed  gases,  with  60%  H2,  40%  N2/CO2,  50  ppm  CO  on 
dry  base,  flow  velocity  =  1  m/s,  temperature  =  82  °C,  pressure  =  3  atm. 


Another  interesting  new  phenomenon  is  that,  over  the  col¬ 
lector  plate  shoulder,  the  chemical  reaction  is  actually  greater 
than  that  over  the  flow  channel.  This  result  is  again  differ¬ 
ent  from  what  was  true  for  a  fuel  cell  using  pure  hydrogen. 
Under  the  collector  plate  shoulder,  both  the  concentrations 
of  CO  and  hydrogen  are  lower.  As  the  mass  diffusivity  of 
CO  is  lower  than  that  of  hydrogen,  over  the  collector  plate 
shoulder,  the  reduced  poisoning  effect  of  CO  is  predominant 
over  the  effect  of  hydrogen  concentration  reduction.  This  phe¬ 
nomenon  is  shown  in  Fig.  7(a)  and  (b).  Fig.  7  shows  anode 
overpotentials  and  power  densities  averaged  over  the  collec¬ 
tor  plate  shoulder  and  over  the  gas  channel,  compared  with 
the  averaged  values  over  the  cell.  This  result  demonstrates 
that  the  shoulder  widths  of  the  collector  plate  at  the  anode 
side  should  be  wider  for  fuel  cells  using  reformate  than  those 
using  pure  hydrogen. 

A  comparison  of  the  results  from  ID  and  3D  model  are 
shown  in  Fig.  8.  In  the  ID  model,  as  the  variation  of  CO 
content  is  neglected  and  an  overestimation  of  CO  poisoning 
effect  results. 


Fig.  8.  Comparison  of  modeling  results  with  ID  and  3D  models,  (a)  current 
density  vs.  overpotential  (b)  polarization  curve.  Base  case  conditions:  anode 
flow  is  fully  saturated  reformed  gases,  with  60%  H2,  40%  N2/CO2,  50  ppm 
CO  on  dry  base,  flow  velocity  =  1  m/s,  temperature  =  82  °C,  pressure  =  3  atm. 
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Fig.  9.  Polarization  curves  at  different  anode  flow  velocities.  Case  conditions:  anode  flow  is  fully  saturated  reformed  gases,  with  60%  H2, 40%  N2/CO2,  50  ppm 
CO  on  dry  base,  temperature  =  82  °C,  pressure  =  3  atm. 


Fig.  9  shows  the  results  of  a  fuel  cell  performance  at  dif¬ 
ferent  anode  flow  velocities.  The  fuel  cell  performance  is 
much  worse  at  higher  anode  flow  velocity.  Generally  speak¬ 
ing  fuel  cell  performance  increases  with  increasing  stoic  ra¬ 
tio,  i.e.  flow  rate,  especially  at  the  cathode.  For  fuel  cell  using 
reformate  as  the  fuel,  higher  flow  rate  at  anode  will  bring  in 
more  CO  as  well  as  more  hydrogen.  Since  under  certain  con¬ 
ditions,  the  effects  of  more  CO  can  overshadow  the  effect 
of  more  hydrogen,  the  resulting  fuel  cell  performance  can 
decrease  with  the  increase  of  anode  flow  rate.  Fig.  9  shows 


Fig.  10.  Polarization  curves  with  varied  anode  flow  velocities.  Case  con¬ 
ditions:  anode  flow  is  fully  saturated  reformed  gases,  with  30%  H2,  70% 
N2/CO2,  50 ppm  CO  on  dry  base,  temperature  =  82  °C,  pressure  =  3  atm. 


such  an  example,  where  the  anode  flow  velocity  increases 
from  1  to  3  m/s  while  all  other  conditions  remain  the  same, 
the  fuel  cell  performance  decreases  considerably.  Therefore, 
under  certain  conditions,  reducing  the  flow  velocity  at  the 
anode  can  be  beneficial.  This  phenomenon  is  similar  to  what 
was  discovered  by  Zhang  et  al.  [21].  However,  there  is  also 
a  limit  to  this  flow  rate  reduction  benefit.  If  the  anode  flow 
rate  is  too  low,  the  limit  current  density  will  be  reduced  as 


Fig.  11.  Power  output  vs.  current  density  with  different  gas  diffuser  porosi¬ 
ties.  Case  conditions:  anode  flow  is  fully  saturated  reformed  gases,  with  60% 
H2, 40%  N2/CO2, 50  ppm  CO  on  dry  base,  flow  velocity  =  1  m/s,  temperature 
=  82  °C,  pressure  =  3  atm. 
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Fig.  12.  Polarization  curves  with  different  width  of  the  collector  plate  shoul¬ 
der  at  anode.  Case  conditions:  anode  flow  is  fully  saturated  reformed  gases, 
with  60%  H2,  40%  N2/CO2,  50  ppm  CO  on  dry  base,  flow  velocity  =  1  m/s, 
temperature  =  82  °C,  pressure  =  3  atm. 

shown  in  Fig.  10.  In  Fig.  10,  the  reformed  gas  contains  only 
30%  hydrogen.  The  anode  flow  velocity  is  3,  1  and  0.5  m/s, 
respectively.  Under  such  conditions,  as  the  hydrogen  supply 
is  already  very  low  and  its  further  reduction  leads  to  reduced 
limit  current  density. 

Beside  the  anode  flow  velocity,  parameters  that  affect  the 
CO  poisoning  include  the  porosity  of  the  anode  gas  diffuser 
and  the  width  of  the  collector  plate  shoulder.  While  the 


cathode  side  benefits  from  narrower  shoulder  and  higher  gas 
diffuser  porosity  due  to  reduced  the  mass  transport  resis¬ 
tance,  same  may  not  be  true  for  the  anode  side.  Fig.  1 1  shows 
the  power  density  of  fuel  cell  with  three  different  effective 
porosities  of  the  anode  gas  diffuser.  When  the  porosity  is 
reduced  from  0.4  to  0.2,  the  anode  overpotential  is  reduced 
and  better  performance  is  achieved,  but  there  is  a  limit  to 
this  reduction.  If  the  porosity  is  too  low,  the  diffusion  rate 
of  hydrogen  to  the  reaction  sites  will  be  very  limited  and  the 
fuel  cell  performance  will  decrease,  as  shown  in  an  extreme 
case  in  Fig.  1 1  when  gas  diffuser  porosity  is  only  0.08. 

Fig.  12  shows  the  modeling  results  with  different  collector 
plate  shoulder  widths  at  the  anode.  For  fuel  cells  using  pure 
hydrogen,  the  wider  the  shoulder,  the  lower  the  performance, 
due  to  the  increased  mass  transfer  resistance.  For  fuel  cells 
using  reformate  with  traces  of  CO,  the  increased  shoulder 
width  will  lead  to  lower  CO  concentration  in  the  catalyst  layer 
and  better  performance.  Similar  to  the  case  with  gas  diffuser 
porosity,  there  is  a  limit  to  this  increase.  If  the  shoulder  is 
too  wide,  the  diffusion  rate  of  hydrogen  to  the  catalyst  layer 
above  the  collector  plate  shoulder  will  be  too  low  and  the  fuel 
cell  performance  will  decrease,  as  shown  in  Fig.  12. 

It  is  obvious  that  the  optimal  design  and  operating  param¬ 
eters  for  a  fuel  cell  using  reformate  could  be  very  different 
from  those  for  fuel  cells  using  pure  hydrogen.  For  a  given 
case,  a  careful  study  should  be  performed  to  select  the  opti¬ 
mal  values. 

Finally,  a  comparison  between  the  modeling  results  and 
the  fuel  cell  performance  data  presented  by  Divisek  et  al. 
[19]  is  plotted  in  Fig.  13,  where  the  computational  results 
compared  well  with  the  experimental  data. 
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Fig.  13.  Comparison  of  modeling  results  with  experimental  data  by  Divisek  et  al.  [19]. 
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4.  Conclusions 

Based  on  our  early  fuel  cell  models,  a  three-dimensional 
mathematical  model  of  PEM  fuel  cell  operated  on  reformate 
is  developed  by  incorporating  the  adsorption  and  oxidation 
kinetics  of  CO  on  platinum  surface.  This  model  is  capable 
of  studying  the  effect  of  CO  poisoning  as  well  as  the  effect 
of  hydrogen  dilution  due  to  the  inert  gases.  In  this  paper, 
presented  are  the  simulation  results  of  PEM  fuel  cell  perfor¬ 
mances  operating  on  reformate  under  various  operating  and 
design  conditions.  The  following  can  be  concluded  based  on 
the  results  of  the  3D  fuel  cell  model: 

•  If  the  anode  stoichiometric  ratio  is  not  too  low,  fuel  cell 
local  current  density  increases  along  the  flow  direction. 
This  is  due  to  fact  that  down  stream  of  the  channel,  the 
reduced  CO  poisoning  effect  is  predominant  over  the  effect 
of  hydrogen  concentration  reduction. 

•  With  enough  hydrogen  supply,  fuel  cell  performance  in¬ 
creases  with  decreasing  anode  flow  rate. 

•  In  the  regular  range,  fuel  cell  performance  increases  with 
the  decrease  of  anode  GDL  porosity.  The  optimal  GDL 
porosity  is  much  lower  for  fuel  cell  using  reformate  than 
using  pure  hydrogen. 

•  There  is  an  optimal  value  for  the  width  of  the  collector  plate 
shoulder.  Usually,  when  contact  resistance  is  neglected,  the 
narrower  the  width  the  better  the  fuel  cell  performance,  but 
for  fuel  cell  using  reformate,  this  is  not  the  case  due  to  the 
opposite  effects  of  CO  poising  and  reduction  in  hydrogen 
concentration. 

The  modeling  results  demonstrated  the  uniqueness  of  3D 
models  over  ID  and  2D  models.  The  results  from  the  3D 
model  have  revealed  new  phenomena  that  cannot  be  predicted 
by  the  previous  ID  or  2D  models.  Please  note  that  the  mod¬ 
eling  results  depend  on  the  design  and  operating  parameters 
and  could  be  different  from  any  specific  set  of  experimental 
results.  Furthermore,  the  core  of  the  3D  model  is  indepen¬ 
dent  of  the  kinetic  models,  and  any  new  kinetic  models  can 
be  easily  incorporated  into  the  3D  model  for  better  results  or 
for  different  catalysts  such  as  Pt/Ru.  The  model  can  also  be 
extended  to  incorporate  effects  of  other  CO  poisoning  reduc¬ 
tion  effects,  such  as  air  bleed,  provided  appropriate  kinetics 
models  are  available. 
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